MAXIMAL REGULARITY WITH TEMPORAL WEIGHTS FOR 
PARABOLIC PROBLEMS WITH INHOMOGENEOUS BOUNDARY 

CONDITIONS 
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Abstract. We develop a maximal regularity approach in temporally weighted L p -spaces for 
vector-valued parabolic initial-boundary value problems with inhomogeneous boundary con- 
ditions, both of static and of relaxation type. Normal ellipticity and conditions of Lopatinskii- 
Shapiro type are the basic structural assumptions. The weighted framework allows to reduce 
the initial regularity and to avoid compatibility conditions at the boundary, and it provides 
an inherent smoothing effect of the solutions. Our main tools are interpolation and trace 
theory for anisotropic Slobodetskii spaces with temporal weights, operator-valued functional 
calculus, as well as localization and perturbation arguments. 



1. Introduction 

In recent years parabolic equations with fully nonlinear boundary conditions have attracted 
a lot of interest since they arise in the analysis of free boundary value problems such as 
the Stefan problem with surface tension, see e.g. [TJ, [11] and |18j . These papers use an L p - 
approach to such problems which yields strong solutions in maximal regularity classes. In this 
framework the boundary conditions are attained in a classical sense up to initial time, and 
not just weakly. This approach is based on linearization and on a sharp L p -regularity theory 
for linear inhomogeneous initial-boundary value problems, as established in [3], [5] and [6] by 
Denk, Hieber, Priiss and Zacher. Besides the usual static boundary conditions, one also has 
to treat dynamical boundary conditions of relaxation type which arise in the context of the 
Stefan problem with surface tension and in related problems. 

However, this approach requires regularities of the initial values (and hence of the nonlinear 
phase spaces) which are stronger than the norms one can control by standard a priori estimates 
for the nonlinear problems. In related situations it is known that one can reduce the required 
initial regularity by means of temporal weights. In the L p -setting, it is natural to work in 

L P ^(J] X):={u:J^X : t^u G L p (J; X)} (1.1) 

endowed with its natural norm, where p G (l,oo), fi G (1/p, 1], T G (0, oo], J := (0,T), and 
t denotes the function t i— > t l ~^u(t) on J. The corresponding weighted Sobolev spaces 
are defined by 

W^(J;X) := {u: J^X : u, u' , G L p>fi ( J; X)} (1.2) 
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for k G N. These spaces and the corresponding anisotropic spaces like ¥, Uifl (J) defined below 
are studied by the authors in detail in |16j . To a large extent they enjoy analogous properties 
as the corresponding unweighted spaces. 

To see the effect of the weight, we consider a generator —A of an analytic semigroup on a 
Banach space X. Then the orbit u(t) = e~ tA UQ belongs to the 'maximal regularity space' 

W^(M+;X)nL PtlM (R + ;D(A)) 

if and only if the initial value belongs to the real interpolation space 

(X,D(A))^ 1/p , p , 

see e.g. Theorem 1.14.5 in |21| . Recall that one often fixes a large p G (l,oo) to treat non- 
linearities. Hence, in the unweighted case p = 1 the resulting initial regularity is close to 
D(A). On the other hand, taking p near 1/p one almost reaches the base space X. Further, 
for Banach spaces X of class HT (see Section 2), Priiss and Simonett have proved in |17j that 
the inhomogeneous evolution equation 

u'{t) + Au{t) = f{t), t>0, u(0) = 0, (1.3) 

has a unique solution in W p ^(M, + ; X) H L P)/i (R_|_; D{A)) for each / G L P:fl (J; X) if and only if 
this fact holds for the unweighted case p, = 1. Since the unweighted case is well understood, 
see e.g. |13j . the L p ^- approach is quite convenient for parabolic problems covered by (|1.3|) . 

Unfortunately, it seems that a sharp regularity theory for inhomogenous boundary value 
problems is not possible within the abstract framework of an evolution equation like (|1.3p . 
Instead one has to restrict to a PDE setting. So we investigate vector-valued linear parabolic 
systems with inhomogeneous boundary conditions, such of static type, i.e., 

dtu + A(t, x, D)u = f(t, x), x G ft, t G J, 

Bj(t, x, D)u = gj(t, x), x G T, t G J, j = l,...,m, (1-4) 
u(0,x) = uq(x), x G ft, 



as well as such of relaxation (or dynamic) type, i.e., 





dtu + A(t, x, D)u 


= f(t,x), 


x g n, 


t G J, 


d t p + B (t,x, 


D)u + C {t, x, Dr)p 


= 9o(t,x), 


x g r, 


£ G J, 


Bj(t,x. 


D)u + Cj(t, x, Dr)p 


= 9j{t,x), 


x g r, 


i G J, 




it(0, x) 


= uo{x), 


x g n, 






p(0,x) 


= Po(x), 


x g r. 





It is assumed that ft C M n is a (inner or outer) domain with compact smooth boundary 
r = (90. In (jl.4p and (|1.5|) the unknown n = u(t, x) takes values in a Banach space i£, and in 
(|1.5p the additional unknown p = p(t,x), which only lives on the boundary V, takes values in 
another Banach space F. Throughout we assume that E and F are of class TiT; for instance E 
and F can be finite dimensional leading to usual parabolic systems. The differential operator 
A is of order 2m, where m G N, and Bj are corresponding boundary operators of order rrij not 
larger than 2m — 1. In (jl.5p the differential operators Cj contain tangential derivatives of any 
order up to kj G No- We assume certain ellipticity and Lopatinskii-Shapiro type conditions 
and impose regularity conditions on the coefficients that are appropriate for the applications 
to quasilinear problems, see e.g. [2], [15]. The details are described in Section 2. 
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We look for strong solutions u of (|1 .4|) . resp. (u, p) of (|1.5p . which satisfy the respective 
equations pointwise almost everywhere. In particular, u shall belong to 

E U)M (J) := W^{J;L p (n;E))nL p ^(j;W^ m (n;E)). 

The space E PiP (J) for p is chosen in accordance to the structure of (|1.5p and to the trace 
theorems established in our paper [16] . see Section 2. In our main results Theorems 12,11 and 
12.21 we show the maximal L Pj/i -regularity for (|1.4p and (|1.5p on J. This means that there are 
data spaces £> s tat(<-0 and T> re \(J) such that (jl.4p and (jl.5p have unique solutions u £ E Mj ^(J) 
and (u,/?) € E MjPj (J) x Ep jAt (J), respectively, if and only if the data satisfies 

(f,9l, -,9m,Uo) G ^stat(^) and (f,go,gi, -,gm,Uo,Po) G £>rel(-0, 

respectively. The data spaces contain the necessary regularities for the data and their com- 
patibility conditions at t = enforced by the static and dynamical boundary equations in 
(|1.4p and (jl.5p . The precise formulations of these spaces is suggested by the space-time trace 
theorems from |16j . In the unweighted case p = 1 and with essentially the same assumptions, 
the maximal L p -regularity for (jl.4p and (jl.5p has been shown by Denk, Hieber &; Priiss [5] 
and Denk, Priiss &: Zacher 0, respectively. 

We note that the problem ([1.5p is more involved in the several respects. Clearly, it contains 
a second variable and a second evolutionary equation. Moreover, the operators Cj can make 
the main parts of the equations at the boundary highly non homogeneous which then leads 
to a rather sophisticated solution space E Pi/i (J) and to a complicated analysis. It further can 
happen that dtp is continuous in t up to t = so that the dynamical equation for p leads to an 
additional compatibility condition on the regularity of i3o (0, •, D)uq + Co(0, •, D-p)po — 9o(0, •)■ 

The main feature of the weighted approach is the flexibility for the regularity of the initial 
values as p varies in 1] . We describe these properties in more detail at the end of Section 2, 
indicating here the basic points. We can solve (|1.4p and (|1.5p with the Besov regularity uo € 
Bp\p 1 ^ p \Q; E) which approaches L p (Q; E) as p tends to I /p. Moreover, if the initial regularity 
is sufficently low we loose the compatibility conditions such as Bj(0, •, D)uq + Cj(0, ■, D-p)po = 
5,(0,-). Since the weight has an influence only at t = 0, our approach yields an inherent 
smoothing effect for the solutions. In particular, for (|1.4p one can control the norm of u(t) 

in Bp^p 1 ^ p \Q;E) by the much lower norm of uq in £p,P 1/p \tt;E). For bounded ft, this 
fact gives the important compactness of the semiflow solving the related nonlinear problems. 
Also for unbounded ft one can thus 'upgrade' the usual a priori estimates in low norms up to 
B^-WfaE) if one is able to handle the involved nonlinearities. See |14j . |15j and also |10| 
in the framework of [17J, as well as [llj. In these papers the weighted approach was used to 
establish convergence to equilibria and the existence of global attractors in high norms. 

In Sections 3 and 4 we prove Theorems 12.11 and 12.21 We first consider model problems with 
homogeneous constant coefficient operators on the full-space M. n and on the half-space Wl 
in Section 3. The full-space case, where boundary conditions are not involved, is treated by 
means of |17] . For the half-space case with boundary conditions we apply the Fourier transform 
with respect to time and space. The solution operators for the resulting ordinary initial value 
problems have been analysed in [5] and [6] for the unweighted case. We now use a recent 
operator-valued Fourier multiplier theorem in the L PjA1 -spaces due to Girardi & Weis [8J and 
several isomorphisms acting on a scale of weighted anisotropic fractional order spaces which 
are investigated in |16] , It is then possible to invert the Fourier transforms and to solve the 
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half-space problem in the required norms. The case of a general domain is then a consequence 
of perturbation and localization arguments, and it is considered in Section 4. 

Finally we discuss several important special cases of (|1.4p and (|1.5p arising as lineariza- 
tions of various quasilinear parabolic problems with nonlinear static or dynamic boundary 
conditions. For instance, the linearization of a reaction-diffusion system with nonlinear Robin 
boundary conditions is of the form (|1.4p for 

E = R N , m=l, A(D) = —A, Bi(x,D) = d u :=vtrnV, 

where A is the Laplacian and v denotes the outer unit normal field on T. The linearization of 
Cahn-Hilliard phase-field models leads to similar problems of order 4 (i.e., m = 2). If we take 

Bi = trn, Ci = -1, 

the static boundary condition for j = 1 in (|1.5p reads u\r = p, which leads to inhomogeneous 
dynamic boundary conditions. Hence p is simply the spatial trace of u in this case. Now one can 
take Cq(x,Dy) = — Ap, the Laplace-Beltrami operator on T, to obtain boundary conditions 
describing surface diffusion, i.e., 

dtu\r + d u u — Ap?x|r = 9o on T, t € J. 

If we choose 

u\r + Ar/> = g\ on T, t e J, 

as the first static boundary condition in (|1.5p . we arrive at the linearization of the Stefan 
problem with surface tension as studied in [7]. Here the graph of p(t,-) is related to the 
unknown boundary at time t. We refer to Section 3 of [6j for further interesting problems that 
may be written in the form (|1.5p . 

Notations. We write a < b for some quantities a, b if there is a generic positive constant 
C with a < Cb. If A is a sectorial operator on a Banach space E, 9 € (0, 1) and q S [1, oo], 
then we set Di{9,p) := [E, D(L)) „ for the real interpolation scale between E and D(L). If 
X, Y are Banach spaces we denote by B(X, Y) the space of bounded linear operators between 
them, with B{X) := B(X,X). 

2. The assumptions and the results 

Throughout we assume that the Banach spaces E, F are of class HT (or, equivalently, are 
UMD spaces). This means that the Hilbert transform is bounded on L2(M;X) which holds, 
e.g., in Hilbert spaces X or if X is a reflexive Lebesgue or (fractional) Sobolev space; see 
Sections III. 4. 3-4. 5 of [1]. We first describe the differential operators in (|1 .4[) and (|1.5p in 
detail. For both problems the operator A is given by 

A(t,x,D)= ^2 a a(t,x)D a , t e J, 

\a\<2m 

where m £ N and D = — iV with the euclidian gradient V = (d Xl , d Xn ) on W 1 . Hence 
the order of A is 2m. The coefficients take values in the bounded linear operators on E, i.e., 
a a (t,x) G B(E). Also for both problems the boundary operators Bj are of the form 

Bj(t,x,D)= bjp(t,x)tt n D p , xeT, tsJ, j = Q,...,m, 

\P\<m, 
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where rrij G {0, ...,2m — 1} is the order of Bj and the coefficients satisfy bo8(t,x) G B(E,F) 
and bjp(t, x) G B(E) for j = 1, m. We note that Bj acts on u by applying first the euclidian 
derivatives and then the spatial trace tr^- We assume that each of these operators is nontrivial, 
i.e., Bj ^ for all j. 

In problem (|l,5p . the boundary conditions of relaxation type involve another set of operators 
Co, ...,C m , which act only on p in the following way. For t G J it is assumed that Cj(t, -,-Dr) is 
a linear map 

C°°(r;F) Li(T;F) 

such that for all j = 0, m, all local coordinates g for T and all p G C°°(T; i* 1 ) it holds 

(c J (t,;D r )p)o g ( x )= ^ cf 7 (t,x) J D^ 1 (po g )(x), xeg-^rnc/), tEJ, 

where [/ C M n is the domain of the chart corresponding to g. Here we have D n _\ = — iV n _i 
with the euclidian gradient V n _i on R n . The order kj G No of Cj is given independently of 
the orders of A and the Bj. The local coefficients cf„ , that may depend on the coordinates g, 
are assumed to satisfy Cq 7 (*,x) G ^(i 71 ) and c| 7 (t,x) G B(F,E) for j = 1, ...,m. It is assumed 
that at least one operator Cj is nontrivial. If an operator Cj is trivial, i.e., Cj = 0, then we set 
kj := — oo as its order. Note that we do not assume that an operator Cj has global coefficients, 
in the sense that there are functions c J7 on T satisfying c? = Cj 7 o g in all coordinates g. In 
contrast to that, the coefficients of the Bj are globally defined on T. The standard examples 
for such an operator Cj are the Laplace-Beltrami operator Ar and a convection term V • Vr, 
where V is a tangential vector field and Vr is the surface gradient on V. Throughout we let 

p G (1, oo) and /U G 1]. 

We look for solutions u of fjl -4|) and (u, p) of (|1.5|) such that 

u G E MiM = J; L P (Q; E)) n ( J; H^ 2m (fi; S)) . 

The weighted vector-valued L Pj/i -spaces and the corresponding Sobolev spaces spaces Wp „ are 
defined in (jl.ip and (|1 .2|) . respectively, and VFp m (0; E) is the valued Sobolev space of order 
2m over fi. For such solutions it the differential equation on the domain fi holds for a.e. (t, x). 
The regularity of u yields 

/ G Eo iM := L P:fl (J; L p (Q; E)). 
From the mapping properties of the temporal trace described in Theorem 4.2 of |16j . we deduce 

u\ t=0 = u e X u ^ := B^-^iQ-E). 

Here Bp 1 (Q;E) denotes the E- valued Besov spaces over fi. We refer to [2], |19| or |24| for 
a definition and the properties of these spaces. Further, Lemma 3.4 of [16J shows that the 
operator maps K u n continuously into 

H l p^ l2m {j-L p {n-E)) nL p ^{J;W2 m - mi (n;E)) (2.1) 

for |^| < rrij < 2m — 1. Due to Theorem 4.6 of [16J, the spatial trace trn maps the space in 
(I2.ip continuously into 

W^{J; L P (T; E)) D L P4l (.J- wf 1 ^ (r ; E)) , j = 0, m, 
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where the number Kj is defined by 

rrij 1 

Kj := 1 , j = 0, ...,m. 

2m 2mp 

Below we assume that Kj ^ 1 — fj, + 1/p for all j = 0, ...,m. The weighted Sobolev spaces 
and the Slobodetskii spaces W£ of fractional order s > are defined by complex and 
real interpolation, respectively. The properties of Wp and Hp are studied in |16j . Moreover, 
Wp(T;E) is the -E-valued Sobolev-Slobodetskii space of order s, where W p = Bp p for s £ No- 
Since the dynamic boundary equation in ([1.5p takes place in F and the static boundary 
equations in (|1.4|) and (|1.5|) take place in £7, these considerations suggest that we should look 
at boundary data 

50 S F 0iM := W£°(j;£ p (T;F)) n^^^^F)), 

g 3 G F it „ := L p (r ; £)) n L p , p ( J; (T; E)) , j = l 



, m. 



For convenience we write 

F p := F , M x ... x F m>M , 5 = ( 5o , ...,5m) G F M , 

and similiarly F p := F iiP x ... x F m , p and <? = (g 1 , ...,g m ) G F M . 

We now determine the regularity of p and po m fll-5p . Assuming sufficient smoothness of 
the coefficients of the operators, we look for a space E PiP for p such that each term in (|1.5p 
involving p acts continuously from E PjP to the space Fj u where the respective equation takes 
place. It can be seen as in Section 2 of [6j that 

=W^(J;L P (T;F)) ni p ,,(J; W l p +2mKo {Y;F)) (2.2) 

satisfies these requirements, where we put 

J := \j G {0, .... m} : A;,- 7^ — 00) L := L -m, + mo, / := max 

1 j=0,...,m J 

It is important to note that 

2m«j + kj = 2rriKQ + j = 0, . . . , m. (2-3) 

We represent E PjP by the points (0,1 + Kq), (I + 2mKo,0), (2mKQ,l) and (kj,Kj), j G J, 
corresponding to the space-time differentiability of the spaces Z% on the right-hand side of 
(|2.2p . The Newton polygon AfV for E p „ is then defined as the convex hull of these points 
together with (0,0). The leading part of AfV is the polygonal part of its boundary connecting 
(0, 1 + Kq) to (I + 2mKo, 0) anti-clockwise. 

Let Zi and Zj be two different spaces on the right-hand side of (|2,2p . It is shown in Propo- 
sition 3.2 of [16J that Z{ H Zj embeds into all spaces whose space-time regularity corresponds 
to the line segment connecting the two points that represent Z{ and Zj in AfV . Consequently, 
the description of E PjP given in (|2.2p contains redundant spaces, in general. We derive a nonre- 
dundant description of E p>p as in the case p = 1 presented in [6] . Here one has to distinguish 
three cases. In each case, a direct application Theorem 4.2 of [16] further yields the temporal 
trace space of p at t = 0, denoted by 



X p ,n '■— tr t= oE P]P . 
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In the same way we also obtain that the temporal derivative dtp has a trace at t = if 
kq > 1 — /j, + 1/p. We denote the resulting trace space by 

Xd tP ,n ■= tr t=0 <9tE P)M if k > 1 - fj, + 

We remark that Theorem 4.2 of [16j means that these trace spaces are given by B pp (T; F) for 
the numbers a > such that (a, k + 1 — [i + belongs to leading part of MV for k = and 
fe = 1, respectively. We can now give the nonredundant description of the spaces M„ „, ^p,^ 
and -Xg 4ftAt . 

Case 1: i = 2m. One has 

Ep,m = < fi0 (J;i P (r;F))nV(j;^ 1+ ^(r;F)) 

since all other spaces in (|2.2p correspond to points on or below the straight line s = 1 + kq — 
r/2m from (0, 1 + kq) to (2m + 2m,Ko, 0) due to (|2.3p . It follows that 

X p ^ = Blf^- 1 ^{Y ] Fl X dt ^ = S M«o-(i-H-i/p))(r ;F ) if k > 1 - /i + 1/p. 

Case 2: i < 2m. One has 

= <+"°(J;L p (r;J0) nL^(j;V^ +2mK °(r;F)) n ^( J; H^ 2mK ° (r ; F)) 

since (1, 2mKo) lies above the line segment s = 1 + kq — r(l + Ko) /(I + 2mno) from (0, 1 + Ko) 
to (2m + 2mKQ, 0) and all points (k^, Ay) are below the line s = 1 + (2m,KQ — r)/l connecting 
(l,2m/«o) an d (0,/ + 2mKo). We then obtain the trace spaces 

= B 2 p ^+i^-m {T . F)) Xdtp ^ = B M«o-(i-m+i/p))(T; F) if K > 1 - /x + 1/p. 
Case 3: I > 2m. Now (1, 2m,Ko) belongs to the interior of MV and it holds 

Ep,m = < K0 (J;i P (r;F))nL w (J;^(r;F)) n n.^^(J;^(r ; F)), 

where J" = {ji, • jq max } C (/max € N, contains those indices j G c7 so that (fej, Kj) belongs 
to the leading part of the Newton polygon, i.e., the points 

P = (0, 1 + kq), Pi = (Aft , Kjl ) , . . . , P 9max = (/c j?max , K j?max ) , F 9max+1 = (Z + 2mK , 0) 

are the vertices of the leading part. It is assumed that J is ordered so that kj q < kj q ^ and 
Kj > Kj for gi < q2- We also define k-\ := 0, := 1 + ko> JTO—i := mo — 2m and 
Z_l := 2m. We further denote the edge in the Newton polygon connecting the points P q and 
P q+1 by 7VP g , q = 0, ...,g max , and set 

:= {j G JU {-1} : (/%■,/%) = P q }, q = 0,...,q mgx , 

J2q+l :={j£jU{-l} : (kj, Kj) G NV q ~), q = 0, g max . 

The temporal trace space of is obtained by Theorem 4.2 of |16) from the spaces corre- 
sponding to Pq and P\. We thus deduce 

X dt ^ = B^ {K0 - {1 -» + ^^ if „ >1- M + 1/P- 

For -Xp,ii one has to distinguish three more cases. 
Case 3(i): If Kj > 1 - \i + 1/p for all j €. J then 

X pi j, = i^+ 2m ( re o-(l-M+Vp)) (P; F). 

Here we apply Theorem 4.2 of [16] to the spaces corresponding to P gmax and P 9l 
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Case 3(ii): Denote by j qi E J be the smallest index with Kj q > 1 — /i+ 1/p, and by j q2 E J 
the largest index with Kj < 1 — fi+ 1/p. Using the spaces corresponding to these indices, we 
conclude that 

% gl +o%- gl -(i-M+vp)) !f :!f 

x™ = Bp# 92 J?1 (r;F). 

Case 3(iii): If Kj < 1 — pt + 1/p for all j € JT, then we employ the spaces related to Po and 
P\ to derive 

fe J1 (K +M-l/p)/(l+/co-/tj 1 ) ,_, 
-&p,fi — &p,p (± , r ). 

For later purpose we finally note that in Case 3 the embedding 

E P ,^W^(J;W;(T;F)), a := J*" , (2.4) 

follows from Proposition 3.2 of [16J by interpolating of the spaces corresponding to Po an d Pi- 
We next consider the assumptions on the coefficients of the operators. For a Banach space 
X of class HT, we write 

F j# (Jxr;l) :=W^{j;L p (T;X))nL p ^(j;Wp mK '(T;X)), j E {0,...,m}. 

The following assumptions shall guarantee that each summand of the operators in (|1 .4[) and 
(|1.5p maps continuously between the relevant spaces described above. In view of the mapping 
properties of the traces and the derivatives, the multiplication with the coefficients has to be 
a continuous map on ¥j jfM (j x T;X) for the relevant X. Moreover, to perform localization 
and perturbation, we require that the top order coefficients of the operators are bounded and 
uniformly continuous. 

(SD) For \ct\ < 2m we have either 2m(fi — 1/p) > 2m — 1 + n/p and a a E Eo,u( J x f2; 13(E)), 
or a a E -Loo (J x £l;B(E)). For \a\ = 2m it holds a a E BUC(J x H;B(E)). If is 
unbounded then for |a| = 2m in addition the limits a a (t, oo) := limi,,.^^ a a (t, x) exist 
uniformly in t E J. 

(SB) Let f = B(E,F) and £ = £(£). For j = 0, ...,m and \(3\ < m, } it holds either b j(5 E 
C T i' 2mT j(7xr;f) with some Tj > Kj, or bjp E Fj i/U (J xT; 5) and Kj > 

(SC) Let J-q = B(F) and T = £>(-F, E?) and let g : V — > T be any coordinates for T. For 
j = 0, ...,m and I7I < kj it holds either c? E C T i ,2mr i (J x r; £) with some Tj > Kj, 
or c| 7 E F J> (J x F; J") and Kj > 1 - fi + 1/p + §=± 

We discuss these assumptions. First, in (SD) one can relax the boundedness assumptions 
for \a\ < 2m, see e.g. |14| . The fact that 

C r 3 ,2 mT] g x p. B ( X ^ . ¥ _^j x p. X j ^ Fj jM (j X T; X) 

for Tj > Kj can be seen using the intrinsic norm for W^{, and Wp mKj given in equation (2.6) 
in |16| and Section 1 of [2], respectively. If Kj > 1 — fx + 1/p + then Theorem 4.2 of |16j 
and Sobolev's embeddings show that 

¥ j;fl (j xT;X) ^ BUC(1 xT;X). 

Using this fact and again the intrinsic norms of Wp/^ and Wp mK:i , we then derive 

¥j ilx (j xT-B(X)) -¥j tlM (J xT;X) F J> (J x 



The assumption Kj > 1 — \i + 1/p + §^ is only valid if p and fj, > l/p are sufficiently large. 
In fact, the assumption is equivalent to p(2m/j> — rrij) > n + 2m. The conditions in (SB) and 
(SC) are not optimal. For all p G (l,oo), one can determine weaker Sobolev regularities for 
the coefficients than the ones given here which still meet the requirements described above, 
see [5], [6] and Section 1.3.4 of |14| . On the other hand, (SB) and (SC) are already sufficient 
for the applications to quasilinear problems, see e.g. [14J and [15J. 

We next state the structural assumptions on the operators, which are the same as in [5] and 
[B]. In the sequel, the subscript ft denotes the principle part of a differential operator, with an 
important exception for the Cj where we put 

C jt := if j i J. 

We thus consider only the principle parts of the operators Cj corresponding to a point on the 
leading part of the Newton polygon for E„ „. First, we assume that A is normally elliptic: 

(E) For all t e J, x eD, and |£| = 1, it holds <r(A$(t,x,t;)) C C+ := {Rez > 0}. If is 
unbounded, then it holds in addition a(A«(t, oo, £)) C C+ for all t G J and |£| = 1. 

We further need conditions of Lopatinskii-Shapiro type. In their formulation, local coordi- 
nates g for the boundary T are called associated to x G T if the corresponding chart (U, (p) 
satisfies 

tp(x) = o, <p'{x) = o u(x) , ip(u nn)ci;, </>(£/ nr) c R^ 1 x {o}, 

where O u m is a fixed orthogonal matrix that rotates the outer normal v(x) of V at x to 
(0, ...,0, — 1) S W 1 . It is easy to see that such a chart (U,(p) always exists. For coordinates g 
associated to x € T, we define the rotated operators A u and Bj by 

A"(t, x, D) :=A{t,x,0^ x) D), B»(t,x,D) := Bj(t, x, O t v{x) D), j = 0,...,m. 

Moreover, we introduce the local operators Cj with respect to g by setting 

Cf(t,x,D n ^):= £ c^g" 1 ^))^, j = 0,...,m, 

\~(\<kj 

where cf., are the local coefficients from the definition of Cj. We continue with the second 
structural assumption concerning (|1.4p . 

(LS s tat) For each fixed t G J and x G T, for each A S and ^ G W 1 ' 1 with |A| + |^'| / and 
each h G -E 1 " 1 the ordinary initial value problem 

\v(y)+A;(t,t',D y )v(y)=0, y > 0, 

B iti(*>£'> A/M^o = hj, j = l,...,m, 

has a unique solution u G Cb([0, oo); E 1 ). 

Here the space Co([0, oo); i?) consists of the continuous E- valued functions on [0, oo) vanishing 
at oo. For the problem (|1.5p with relaxation type boundary conditions we need two assumptions 
of Lopatinskii-Shapiro type in the Cases 2 and 3. First, in each case we require a natural 
analogue of (LS stat ). 
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(LS re i) For each fixed x G T, choose coordinates g associated to x. Then for every t £ J. 

A G C+ and £' G W 1 ^ 1 with |A| + ^ 0, h £ F and € E, j = l,...,m, the 
ordinary initial value problem 

(\ + A u t (t,x,Z',D y ))v(y) = 0, y>0, 

^(M^A/Mi^o + {\ + C^(t,x,^'))a = h , 

Bj t (t,x,£',Dy)v\ y= o +C^(t,x,£')a = hj, j = l,...,m, 

has a unique solution (v, a) £ Co([0, oo); E) x F. 
In the Cases 2 and 3, the following additional 'asymptotic' conditions are required, respectively. 
(LS^) Let I < 2m. For each fixed x £ T, choose coordinates g associated to x. Then for every 

t G J, ho £ F, hj EE, 3 = 1, ...,m, and each A G CjT, f G R n_1 with |A| + 7^ 0, 

the ordinary initial value problem 

(X + A^(t,x^',D y ))v{y) = 0, y>0, 

Bj$(t,x,€',D y )v\y=o = hj, j = l,...,m, 

and for all A G C+ and = 1 the problem 

A\{t,x,t',D y )v{y) = 0, y>0, 

Bo^(t,x,^,Dy)v\ y= o+ {X + C^(t,x,^'))a = h , 

Bj$(t,x,£',Dy)v\ y= o+C^(t,x,£')a = hj, j = l,...,m, 

has a unique solution (v,a) G Co([0, 00); E) x F, respectively. 
(LS^) Let I > 2m. For each fixed x G T, choose coordinates g associated to x. Then for each 
t G J, ho G F, foj £ E, j = 1, ...,m, and each A G C + and £ G M n_1 \{0}, the ordinary 
initial value problem 

(A + ^(t,x,^,D y ))t;(y) = 0, y > 0, 

^. -Ot/)^!!/^ + ^ Jw+i C l(*> I » O ff = h v j = 0,...,m, 

and further for all A G C+\{0}. = 1 and 9=1, ...,2(7 max , the problem 

(A + ^^^O,^))^^) = 0, y>0, 

Boj t (<,a;,0,Dj / )u[ y= o + 5-i )> 7 9 Acr + (yo )> 7 9 Cf )t (i ) a; ) £ / )o- = ^0, 

BJjt(t,x, 0,^)^1^=0 + j^Cljj^x,^')^ = hj, j = l,...,m, 

has a unique solution (v,a) G Co([0,oo);F) x F, respectively. Here 5^ y = 1 if j G Jq 
and <5j ; j- = otherwise. 

In [5] and [6] it is shown that these conditions are necessary for maximal L p -regularity of (jl.5p 
on finite intervals. In Section 3 of [6| they are verified for a variety of concrete problems from 
the applications, see also [14] and |15] , If F and F are finite dimensional, it suffices to consider 
ho = hj = in the above conditions. 

We can now state our maximal L Pi/J -regularity results. We start with the one for (|1.4p . 

Theorem 2.1. Let E be a Banach space of class "HT, p G (1, 00) and p G (l/p,i\. Let 
J = (0, T) be a finite interval, and let £1 C M. n be a domain with compact smooth boundary 

10 



r = dil. Assume that (E). (LS s t a t)> (SD) and (SB) hold true and that Kj 7^ 1 — fi + 1/p for 
j = l,...,m. Then the problem 

dtu + A(t,x, D)u = f(t,x), x G ft, t G J, 

Bj(t, x, D)u = gj(t, x), x G T, t & J, j = l,...,m, 
u(0,x) = uo(x), x G 0, 

has a unique solution u = C sta ,t(f,g,uo) G if and ouiy if (f,g,u ) G P s tat, where 

£>stat := {(f,9, uo) G E , M x F p x X UiP : for j = 1, m it holds 

Bj(0, ■, D)uq = gj(0, •) on F if Kj > 1 — \i — 

The corresponding solution operator £ s tat : ^stat E UjAt is continuous. If £ s tat is restricted to 

^stat : = {(/, <7,^o) e ^stat : ffj|t=0 = if Kj > 1 - ii - j = 1, ...,77t}, 

for any given To > the operator norm of the restriction is uniformly bounded for T G (0, To]. 

In the situation of the theorem, it is clear that if the coefficients 

(-i) |a| a Q , \a\<2m, (-i) |/3| 6 j/3 , \/3\ < m j} j = l,...,m, 

and the data are real-valued, then also the solution u is real-valued. We next state the maximal 
regularity result for (|1.5p . 

Theorem 2.2. Let E and F be Banach spaces of class HT, p G (1, oo) and p G (1/p, 1]. Let 
J = (0, T) be a hnite interval, and let £1 C W 1 be a domain with compact smooth boundary 
r = d£l. Assume that (E), (LS re l) ; (SD), (SB) and (SC) are vaiid and that, in addition, if 
I < 2m condition (LS^) holds and if I > 2m condition (LS^) holds. Assume further that 
Kj 1 — p, + 1 Jp for j = 0, m. Then the problem 





dtu + A(t, x, D)u 


= f(t,x), 


x G 0, 


t G J, 


d t p + B (t,x, 


D)u + C (t,x,D r )p 


= 9o(t,x), 


x G r, 


t G J, 


Bj(t,x, 


D)u + Cj(t, x, Dr)p 


= 9j(t,x), 


x g r, 


t G J, 




u(0, x) 


= u (x), 


x G O, 






p(0,x) 


= Po(x), 


x g r, 





has a unique solution (u, p) G E UiP x Ep )At if and only if (f,g,uo,po) G T> r el> where 

A-el := {{f,g,u ,p ) G E 0j(U xFpX X U ^ L x X p>p : for j = 1, m it hoids 

B^O, •, £>)«o + Cj(0, ■, D r )p = gj {0, •) on T if Kj > 1 - p. + 1/p; 
50(0, •) — $o(0, ')D)uq — Co(0, •, Dy)pq G Xd tP; n if k q >1- p + l/p). 

The corresponding solution operator C re \ : T> re \ — > E UiP x E p>p is continuous. If C re \ is restricted 
to 

V m\ '■= {(f,9,u ,p ) G V ml : pj| t=0 = if Kj > 1 - p - 1/p, j = 0,...,m}, 

for any given Tq > the operator norm of the restriction is uniformly bounded for T G (0, To]. 
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A similiar statement as above holds for real- valued solutions. In the theorems, the spaces 
Pstat and V Te \ are considered as closed subspaces of Eo,/i E^ x X u ^ and Eo,/i x E^ x X u ^ x X p ^, 
respectively. They contain the compatibility conditions of the boundary inhomogeneities and 
the initial values at t = 0, which are necessary for the solvability of (|1.4|) and (|1.5p . respectively. 

One needs the spaces 2?g tat and Pj? el with vanishing initial values since the resulting uniform 
estimates are crucial for fixed point arguments arising in the context of quasilinear problems. 
They are considered as closed subspaces of Eo : « x q¥^ x X u „ and Eo,^ x qF^ x X u „ x X p n, 
respectively, where oE^ and qF^ are defined as follows. For a Banach space X of class HT and 
s = [s] + s* with [s] £ No, s* G [0, 1) we set 

W;jJ;X) := { W^(J;X), W p f 1 (J;X)) s ^ p , 

where W^(J; X) := {u G W^(J;X) : n(0), n^" 1 ) (0) = 0} is considered closed 
subspace of W p> JJ; X), and then 

oE iiM := W^(J; L p (r ; E)) D (r ; £)) , j = 1, .., m. 

Analogously we define the spaces oEo,^, oE^, oE^, oE M ,^ and oEp,^. It is shown in Proposi- 
tion 2.10 of [16] that W^ = if 6 < s < 1 - p + 1/p and 

0% S ,, = {«6^ : « (0 (0)=0, l€{0,...,k}} 

if — /x + l/p < s < k + 2 — p + 1/p for 6 No, with equivalent norms, respectively. In 

other words, the trace at t = of a derivative of u G oW p vanishes if it exists. 

Compared to the unweighted approach, the maximal L PiAt -regularity approach has the fol- 
lowing advantages, where we restrict to the setting of f)l -4[> . Analogous statements are valid 
for (fl~5l) . 

• Flexible initial regularity: We obtain solutions for initial values in B pp (p,;E), where 
s G (0,2m(l - 1/p)}. 

• Inherent smoothing effect: Away from the initial time, t G (0, T), the solutions belong 
to 

W} (r, T; L p (n- E)) D L p {r, T; W^^; E)) ^ C(7; B^ 1 '^ (Q; E)) , 

• Control solutions in a strong norm at a later time by a weaker norm at an earlier time 
and the data: For s = 2m(p — 1/p) G (0, 2m(l — 1/p)] it holds 

• Avoid compatibility conditions: Given p G (l,oo), if p is sufficiently close to 1/p, we 
have Kj < 1 — p + l/p for all j and thus obtain a unique solution u G E U) „ for arbitrary 
data in E 0)M x x X u ^. 

The rest of the paper is concerned with the proofs of the Theorems 12.11 and 12.21 

3. The model problems 

We first consider the full-space case Vt = W 1 without boundary conditions and assume that 
the coefficients of the differential operator 

A{D)= Y, a ^D a 

\a\=2m 
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are constant, i.e., a a G B(E) are independent of (t,x). Observe that there are no lower order 
terms and that A is homogeneous of degree 2m. We have the following maximal L Pi/J -regularity 
result for A on the half-line. 

Lemma 3.1. Let E be a Banach space of class HT, p G (l,oo), fi G (1/p, 1], and assume 
that the constant coefficient operator A satishes (E). Then there is a unique solution u = 
S F {f,u ) eE v (l + xK") of 

u + d t u + A(D)u = f(t,x), x G R n , t>0, 

u(0,x) =u (x), xeR n , (3.1) 

if and only if 

f € Eo j(U (M+ x R n ) and u eX U)fl (R n ). 

Proof. Lemma 4.2 of [5] shows that the realization of the operator 1 + A on L p (R n ;E) with 
domain D(l + A) = W p m (R n ; E) admits maximal L p -regularity on the half-line. Since 

X^(R») = B^- l M{R n -E) = {L p (R n ;E),W^ m (R n ;E)) fi _ 1/p ^, 

the assertion follows from Theorem 3.2 of [17J. ■ 

The model problems for (|1.4p and (|1.5p on the half-space involve boundary conditions 
and thus require a much greater effort. To construct the solution, one uses an operator- 
valued Fourier multiplier theorem in L p n. For Banach spaces X, Y and a symbol m 6 
£l,k>c(^; B(X, Y")) one introduces an operator T m by setting 

T m f := T^mTf, f G T^C^(R; X), 

where J- denotes the Fourier transform on the real line. We can restrict T m to functions on 
K + . Observe that T m is densely defined on L p ^(R + ; X). We also use the analogous definition 
on the space L p (R n ; X). The next result is due to Girardi & Weis [8]. 

Theorem 3.2. Let p G (l,oo), p, G (1/p, 1], and let X,Y be Banach spaces of class HT ■ 
Assume that m G C 1 (M\{0}; B(X, Y)) satishes K({m{\), Am' (A) : A / 0}) < oo. Then 
T m eB(L p ^(R + ;X),L P)fM (R+;Y)). 

Here, the 7£-bound of a family T C B(X,Y) is denoted by 1Z(T). For a definition and 
properties of 7^,-boundedness we refer to [3] or [13] . Under more restrictive assumptions on the 
symbol m we can give a short proof a multiplier theorem in L p „, employing a result of Kree 
[12J (which is also used in the proof in [8j). 

Proposition 3.3. In addition to the assumptions of Theorem \3.2[ suppose that m satishes 

m G C 2 {R\{0};B(X,Y)), \m"{X)\ B{xy) < \X\- 2 for A + 0. 
Then T m G B(L p ^(R + ; X), L p ^(R + ;Y)) . 

Proof. The operator- valued multiplier theorem for the unweighted case \i = 1 shows that T m 
extends to a bounded operator from L p (R + ; X) to L p (R + ;Y); see Theorem 3.4 of |22j . More- 
over, following the lines of the proof of Lemma VI. 4. 4. 2 of |20j, the assumptions on m imply 
that T m may be represented as a convolution operator with a kernel k G C(M\{0}; B(X, V)) 
satisfying \k(t)\^^x,Y) ~ 1^1 ■ ^ now f° nows from Theoreme 2 of |12j that T m is also bounded 
from L p ^(R + ;X) to L P4l (R + ;Y), for all fi G (l/p,l]. ■ 
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We next treat the half-space model problem corresponding to (|1.5p . where we proceed simi- 
larly as in Section 4 of [6]. On Q = R™ with boundary T = R n_1 we consider the homogeneous 
differential operator 

A(D)= a a D a 

\a\=2m 

and the homogeneous boundary operators 

B j {D) = bjptronDP, Cj{D^ x ) = £ c h Dl_ x j = 0,...,m. 

\/3\=mj h\=kj 

The coefficients of the operators 

a a ,b jP eB{E), Cjl eB{F,E), j = l,...,m, b pEB{E,F), c 07 G B(F) 

are assumed to be independent of t and x. If nothing else is indicated, now all spaces have to 
be understood over R + x R™ and over R + x R n_1 . 

Lemma 3.4. Let E and F be Banach spaces of class TiT, p G (l,oo), and p G 1]. 
We assume that (E) and (LS re i) are vaiid and that condition (LS^) holds if I < 2m and 
condition (LS+ ) holds if I > 2m. Let (f,g,Uo,po) G T> re \. Then there is a unique solution 

(u, p) G E Ujfl x E PiAt of 



u + d t u + A(D)u = f{t,x), 


x G I 


^, 


t > 0, 


p + d t p + B (D)u + C (-D n -l)p = 9o(t, x), 


x G I 


l n ~\ 


t > 0, 


Bj(D)u + Cj(-D„-i)p = gj(t, x), 


x G I 


l n ~\ 


i>0, j = l,...,m, (3.2) 


u(Q,x) = u (x), 


x G I 


11, 




p(0,x) = po(x), 


x G I 






Proof. (I) We first show uniqueness for (|3.2I). We 


use the 


space 


Z := L p (R";^)xW p s (R n - 1 ;F) 



with s = 2mKo in the Cases 1 and 2 as well as s = kj 1 Ko/(l + k.q — Kj x ) in Case 3. On Z, we 
introduce the operator A defined by 

A(u, p) := ((1 + A)u, B u + (1 + C )p) , (u, p) G D{A), 

with domain 

D(A) := {(u,p) G W p 2m (R" ;J B) x Wp +2mK ° (R n_1 ; F) : 

BjU + (1 + Cj)p = 0, j = l,...,m; B « + C p G W^R^ 1 ; F)}. 

By (the proof of) Theorem 2.2 of [6], A generates an analytic Co-semigroup on Z. Due to (|2.2p 
and (12. 4p . the space E u „ x E„„ embeds into 

G :=E n , M (R + ) x (^ )M (R + ;W p s (R n - 1 ;F))nL Pj(U (R + ;^ +2mK0 (R ri - 1 ;F))). 

Let « £ G be a solution of (J372]) with tto = 0, po = 0, / = and go = ■ ■ ■ = g m = 0. Since 
L P u(J',Z) Li(J',Z), it follows that u is a mild solution of the inhomogeneous evolution 
equation for A on Z with trivial data, and thus u = 0. 

(II) The rest of the proof is concerned with the existence of solutions of (|3.2[) . We write 
x = (x',y) G R+ with x' G R n_1 and y > 0, as well as T x < and F( for the partial Fourier 
transform with respect to x' and t G R, with covariable £' G M n_1 and G R, respectively. In 
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order to apply Tt, we extend a function with compact support in R + by to R. In the same 
way as in Section 4.1 of [6] one can see that it sufficies to consider the case 

/ = 0, g = (g , ...,g m ) G oF M , u = 0, p = 0. 

(See Lemma 3.2.2 and Proposition 3.2.3 of [II].) Moreover we first assume that 

geV:= C C °°(R + x R™- 1 ; F X E m ). 

It can be seen as in Lemma 1.3.14 of |14) that T> is dense in qF^. For such data the problem 
(|3.2p was solved in the proof of Theorem 2.1 of In the following we estimate the norm of 
the solution (u, p) in the weighted solution space oK Uifl x o^p,n by the norm of g in qF^. For 
this estimate, we have to derive an appropriate representation of (u,p). We apply T x iTt to 
(|3.2p and arrive for any 9 G R and £' G R n_1 at the ordinary initial value problem 

(l + iO)v + A(g,D v )v = 0, y>0, 

(l + i0)a + Bo(e,D y )v\ y=o + C o (e)a=(F x ,F t g o )(9,g), (3.3) 

B 3 (t',D y )v\ y=Q + C j {0<?= fa?t9j)M), j = l,...,m. 

In Section 4.3 of [6| it is shown that (j3.3|) possesses for all 6 and £' a unique solution 
(v(9,^', -),a(9,£')) which may be represented as follows. We define the symbols 

tf := (l + i0+|£'| 2 ™)V2™ 6:= M C:= 4, a - 1 + ^ 



and the so-called boundary symbol s(0,£') by 

s(9, £') := 1 + i0 + in the Cases 1 and 2, 

s (0,f) --i + W+Y, \t'\ k i-& m °- m i in Case 3. 

Then it holds 

v(9,C',y) = first component of P s (bC, a)M°{b, £<?) [p-^Wtgjtf, O) i=0 ,.„,m' 

= K^cO"^ mo M^6,c^)(^- mj ^^(^0), = o,..., m - 

Here we have used holomorphic functions 

A) : C"- 1 xC4 B(E 2m ), P s : C™" 1 xC^ B(E 2m ), 
M„ D :D t xD ( xE4 B(F x £ m , £ 2m ), M° : £> b X £> c X S -> B(F x E m ,F), 
where Z?b C C and C C n_1 \ {0} are bounded open sets satisfying 

(B 1/2 (l/2)) 1 /2- c Dbf {( e R n-i . | C | = 1} c Do 

and S = S^ = jz G C\{0} : |argz| < 0} is a sector with G (^, 7r). The spectrum of iAo(bC, a) 
has a gap at the imaginary axis, and P s (bC,a) is the spectral projection corresponding to the 
stable part of the spectrum. The maps M® and M® have the crucial property that 

{Wl^D^M^leWr 1 ,^) ■■ a'e{0,l} n -\ CVO, b G D b , d G S} (3.4) 
is an 7^-bounded set of operators in B(F x E m ,E 2m ), and that 

{le'l^'l^'M^b^'lCr 1 ^) : a'G{0,l} n -\ eVO, 6 G D b , •& G S} (3.5) 
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is an 7^-bounded set of operators in B(F x E m , F). For the solvability and the representation 
of the solution of (|3.3p in [6] only the condition (LS re i) is needed. In the Cases 2 and 3 the 
asymptotic Lopatinskii-Shapiro conditions (LS^,) and (LS+ ) are required to show the TZ- 
boundedness of the sets in (|3.4p and (|3.5p . because of the unboundedness of In Case 1 the 
symbols M® and M® do not depend on i?, so that in this case additional conditions are not 
necessary. 

Since T x iTig belongs for g G T> to the Schwartz class and all derivatives of the terms 
involved in the representation of the solution grow at most polynomially, we can apply the 
inverse Fourier transforms and obtain that 

u= first component of F7 1 e^ b ^ P S {K, a)M°(6,C,tf) [fi^'Wig,) i=0 ,..., w > 

is the unique solution of (|3,2p with f = 0, uq = 0, po = and g G D. 

(Ill) We derive another representation of the solution by identifying the Fourier multipliers 
with operators. For a function h G S(W l ^ 1 ;E 2m ) and fixed (x',y) G WL we calculate 

(Tfj» A °yp s h)(x') =(T-ij» A °(y+®p s e-y#h)(x% =0 (3.6) 

/•oo 

= ~ / dy{T-U* Ao{y+ ®P s e- y »h)(x')&y 
Jo 

= jH (Ffe^(y^P s ^^^e-y%)(x')dy 

[•CO 1 ; A 

= l (^ 1 e i ^^P s y^)*((W i - 1 / i )(,y))(,')d^ 

neglecting the arguments of Aq and P s . Here the operator Lq is defined by 

L e := 1 + W + (-A n _!) m = T-}d 2m T xl , 

where the last equality holds, e.g., on Schwartz functions. We observe that for a bounded 
holomorphic scalar function pna sector S r with r G (0,tt) the operator <p(— A n _i) defined 
via the 7^°°-calculus for — A n _i on L p (]R n ~ 1 ; E) coincides with the Fourier multiplier J' x ~/ li p(\ • 
| 2 )Jv, see Example 10.2 of |13j . Moreover, the 'H 00 -calculus extends the usual Dunford type 
calculus for sectorial operators, see Remark 9.9 of [13J. Therefore, the extension operator £q, 
which corresponds to y i— > e -2 ^, is given by 

(£ e f)(x', y) := e-v Ll * 2m f(x'), x' G R n ~\ y > 0, 

for / G L p (M n ~ 1 ; E). We also obtain the equality 

T-}$ 2m e-®h = L e S e T-}h, h G 5(M n_1 ; E m ), 

which we have used in the last line of (|3.6p . For 6 G M and / G ; £' 2 ' m ) we thus define 

the operator T{9) by 

(T(9)f)(x',y) := first component of J (^-^^0 Wp s __0) „ /(., y) ( x ') d y. 

The proofs of Lemmas 4.3 and 4.4 in [5] show that T G C^R; £(L p (R^; E 2m ), W 2m (R%) E))) 
and that 

\D a T(9), 9^-D a T{6) : # G R, lal < 2m) (3.7) 
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is an 7~£-bounded set of operators in £?(L p (R™ ; E 2m ), L p (R™ ; E)^ . Further, as above one can 
see that r d~ m iT x i = T x iL Q m ^ 2m on Schwartz functions, for j = 0, ...,m. This fact leads to 

The Dunford type calculus for sectorial operators yields for 9 € R and y > the representation 

L e e-y L l /2m = ±f s Ze-v* 1/2m (z-L e yUz, 

where S = (oo, <5]e l37r / 2 U Je 1 ! 371 "/ 2 ' -37 ^ 2 ! U [5, oo)e _l37r/ ' 2 for some sufficiently small 5 > 0. Hence 

/IN T ^ /-^ 7n 

for each y > the B(L p (R ri , i?) J -valued function i— > Lge~ y « is smooth and all of its 
derivatives are bounded. So we can apply the inverse Fourier transform with respect to i and 
obtain that 

on Schwartz functions, where L := l + dt + (— A n _i) m and £ := e~' Ll/2m . Here, for X € {i^i 71 } 
we consider L as an operator on L P)At (R + ; L p (W l ~ 1 ; X)) with the domain 

D(L) = D(d t ) + D((-A n . 1 ) m ) = W^(R+;L p (R n - 1 ;X)) nL p ,,(M + ;^ 2m (R n " 1 ;^))- 
In Lemma 3.1 of [16J we have established that L is invertible and sectorial with angle not larger 
than 7r/2. Similarly one can treat fractional powers and derive L m ^ 2m = T t L- m ^ 2m Tf 1 . 
We arrive at 

u = C u g := [Tt l T{9)T t ) LS {JT^Kfo (^ /2m #) i= o,...,m> 
Analogous arguments show that the second component p can be represented by 

p = C p g := S - l L m ^ 2m {T^TfM^b, £ V)F x >T t ) (L~ m J / 2 ™ 9j ) m , 
with the operator 

S ■= 1 + d t + (-A n _i) //2 in the Cases 1 and 2, 
S:=l + d t + J2{-^n-l) kj/2 L^ mo - m ^ 2m in Case 3. 

Using the properties of L proved in Lemma 3.1 of |16j . it can be shown as in Section 4.2 of 
[6] that S is an isomorphism between o^p,p an d 0^0,^- Because T> is a dense subset of oF^, it 
now remains to prove the estimate 

\£ugk u ,„ + \£ P gk P ,K < I#| f m , (3.8) 

If (|3.8p has been verified then the solution operator C := (C u ,C p ) extends continuously to an 



operator from qF^ to qE u ^ x oEp,u, and this extension yields the solution of (|3.2p . 
(IV) Lemma 3.1 of [16] says that for s G (0, 1] we have 

D L (s,p) = %%(M+;L p (]R n - 1 ;X)) nL PiM (M + ;iy p 2ms (M™- 1 ;X)). 

Therefore, for j = 1, m the operator j J - m il 2ra maps the space o^ju = Dl(kj,p) continu- 
ously into 

Y E := D L (1 - l/2mp,p) = W^ 2mp {R + ; L^R"" 1 ; E)) n L Pj(U (R + ; W J 2m ~ 1/ ' p (M n ~ 1 ; E)) . 
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The same arguments yield that L m o/ 2m maps o^o,^ continuously into o^F> which is defined 
as oYe with E replaced by F. We next prove that the operator 

on V with the symbol M° : D b x D c x S ->■ B(F x E m , E 2m x F) given by 

M°(6,C,tf) := (M u (6,C,^),M p (6,C^)), 

extends continuously to an element of B^Yp x oY^^Y 2 ^ x o^f). To this end, we consider 
the approximating operators 

M°' £ ■.= ^ 1 T- l 1 M Q {b,C^){l+$r £ J : x>J r u £6(0,1). 
Observe that Ai°' £ (l + L 1//2m ) e = .M on V. Cauchy's formula yields the representation 

M °' 6 = ~T* [ [_ ^r 1 ^ 1 ^°(^C,^)(i + $r £ (b -b)- 1 ^ - ay^FtdbS, 

with 3$ = ( — oo,0]e _1< ^* U [0, oo)e 1 ^* for some £ (ir/Am, (j)), and where is a closed curve 
in surrounding (S 1 / 2 (l/2)) 1//2m . Since £ = 4^ is independent of 0, we may rewrite the 
above equality as 

_M<^ = _J_ j j FfM°(b,(,d)T x ,(l + d)~ £ (b - B)- l (d - i}l 2m Y l d&dtf, 

where B := (— A n _i) 1 / 2 L _1 / 2m corresponds to the symbol 6 = ^U. The realization of 5 
on Lp ]At (R + ; L p (R n_1 ; is 1 )) is a bounded operator, and its spectrum is contained in the set 
(-B]/ 2 (l/2)) 1 ' 2m . This can be seen using the joint functional calculus for dt and ( — A„_i) m on 
L PiM (]R + ;Lp(]R n - 1 ; J E)), see Theorem 4.5 of 0. 

Due to the 7£-boundedness of the sets (|3.4p and (|3.5p , the operator- valued Fourier- multiplier 
theorem in R n_1 (Theorem 3.25 of [3], see also Theorem 4.13 of [13]) and real interpolation 
imply that that the operators 

M 1 ^, $) := F~, l M°(b, ■,'&)T X >, b € D b , G E, 

extend continuously to elements of S(W^(R n-1 ; .F x £ m ), W^(R n_1 ; £ 2m x F)), s > 0, with 
uniformly bounded operators norms. Since M° is holomorphic, also M 1 is holomorphic in its 
arguments. By canonical pointwise extension we thus obtain that 

M 1 : D b x £ -> £( Y F x Yg, Y| m x Y F ) 

is bounded and holomorphic. Using L as an isomorphism D(L) -> L PlM (R+; L p (R n_1 ; F)) that 
commutes with F, we see that the spectrum of the realization of B on D(L) is also contained 
in (F 1 / 2 (l/2)) 1 / 2m . By interpolation, the same holds on qYp x o^^. Hence, we may rewrite 
M°' £ as 

M°> £ = -^ 1 [ j M 1 (b,d)(l + dy £ (b-B)-\d-L 1 / 2m )- l dbdd, 
where the curve integrals are now defined in B(qYf x o^^oY 2 ^ x oYp). We thus obtain 

M°' £ = — / M 2 (tf)(l + tf)" e (tf -L 1 / 2 ™)^?? 
2vri 7s,, 
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for a bounded holomorphic map 

M 2 : E — )■ 0( O Y F x Y^, Y| m x Y F ). 

Since the realization of L 1//2m on L PjA1 (R + ; L p (R n_1 ; is sectorial with angle not larger than 
7r/4m, it follows from Corollary 1 of [3] that L 1 / 2 " 1 admits a bounded operator- valued Ti 00 - 
calculus with %°°-angle not larger than ir /Am on the real interpolation spaces qY 1 ^ and oY F , 
respectively. From this fact and the boundedness of M 2 on E we infer 

I^°' £ |b(oY f x Y™ Yf"x Y F ) ~ S U P |Af 2 ('!?)(l + ^)~ £ |e(oY F xoY™,oY|™x Y F ) < C> ( 3 - 9 ) 

where C does not depend on e € (0, 1). Due to Proposition 2.2 of [3], for /i € D{L 2 ) the map 
e i— >■ (1 + L 1//2m ) e /i is continuous with values in Dl(1 — l/2mp,p). Together with f|3 . 9[> . this 
fact yields 

l- A ^ 0/l l Y 2 ™x Y F ^ limsup |^ ' £ |b( oYf xoY™,oY 2 ™xoY f ) 1(1 + ^ 1/2m ) e/l loY F xoY™ < N Y F x Y£ • 

£->0 

Since D(L 2 ) is dense in -D_l(1 — l/2mp,p), we obtain that .M extends to an element of 
B( Y F x Y™, Y| m x Y F ), as asserted. 

(V) Now we can show the required estimate for C u , i.e., 

|£u5k iM < \g\ v M , #e£>. (3.10) 

The extension operator £ = e ~- Ll/2m maps continuously 

D L (1 - l/2mp,p) = D L1/2m (2m - l/p,p) -> L P (R+;D(L)), 

and L maps the space L p (R+;D(L)) continuously into 

L p (R + ; L PiM (R + ; L p (R n = Lp iAt (R+; L p (R™ ; E 1 )) . 

Of course, here E may be replaced by F. Thus L £ maps continuously 

Y| m x Y F -»• L P4l (R + ; L p (R^; E 2m x F)) . 

Theorem 13.21 and the 7£-boundedness of (|3.7p imply that J r f 1 T(-)Ft extends to a continuous 
operator 

L p ^{R + ;L p (R n + ;E 2m )) -+ L PiAt (R + ; W 2m (R n + ; E)) . 

Alternatively, this fact follows from Proposition 13.31 since one can show that the operator 
family 

{0 2 ^D a T{9) : \a\ < 2m} 

is bounded in S(L p (R"; E 2m ), L p (R™ ; E)) arguing as in the proof of Lemma 4.4 of [5]. The 
equation for u shows that its E UiA1 -norm can be controlled by its Lp jA1 (R + ; W 2m (R™ ; £'))-norm. 
So we have established (|3.10p . We finally consider the required estimate for C p . As above we 
obtain that L m o/ 2m maps continuously 

Y F = D L (1 - l/p,p) D L { Ko ,p) = F 0)M . 

Since 5 _1 is an isomorphism from qWq^ to o^p,^) this gives the estimate for C p . ■ 

The analogous half-space result for (|1.4p reads as follows. 

19 



Lemma 3.5. Let E be a Banach space of class HT, p G (1, oo), \i G (1/p, 1], and assume that 
(E) and (LS) are vaiid. Then for (f,g,uo) G 2? s tat there is a unique solution u G E U) „ of 



u + d t u + A{D)u = 


= f(t,x), 


x G R™ , 


t > 0, 






Bj(D)u = 


= 9j{t,x), 


s g m n -\ 


t > 0, 


j = l,...,m, 


(3.11) 


u(Q, x) = 


= u {x), 


x G R™ . 






■ 



We refrain from giving a detailed proof of this result, which is similar to the one of Lemma [3.4l 
and also less sophisticated. (See Section 2 of |14| for the details.) Again we may restrict to 
the case / = 0, g G o¥^ and uq = 0. Applying the partial Fourier transforms with respect to t 
and x' to (|3. we arrive at an ordinary initial value problem, whose solution operator is for 
regular data (gi,...,g m ) given by 

due to Lemma 4.3 of [5j. Here T has the same properties as T and L, 8 are given as in the 
proof of Lemma 13.41 The arguments given in the Steps IV and V of the proof above yield that 
C G B(q¥^, oIE Ui ^), which implies the solvability of (|3,lip as asserted. 

4. The general problem on a domain 

Theorems 12.11 and 12.21 are now a consequence of the above results for the model problems 
and a perturbation and localization procedure, analogous to the one in e.g. Section 4.5 of [6]. 
We only sketch the proof below since the full procedure is rather lengthy and tedious. The 
arguments are worked out in great detail in Sections 2.3, 2.4, 3.2.2 and 3.3 of |14j . Moreover, 
we concentrate on (|1.5|) since the proof for (jl.4p is similar and a bit simpler. 

Proof of Theorem 12.21 (I) Let us first consider the necessary conditions on the data. 
The considerations in Section 2 and the assumptions (SD), (SB) and (SC) yield that A G 
B(E U)M ,E 0)At ) and Bj G B(E u ^,¥j^), Cj G #(E PiAt , Fj )At ) for j = 0, ...,m. Moreover, we have 

Wp^(j; L P (T; E)) ^ BUC(1; L P (T; E)) if Kj>l-fi+ 1/p 

for j = 1, ...,m, due to Proposition 2.10 of |16| . Thus in this case the j-th. boundary equation 
in (|1.5p must hold up to t = by continuity, which explains the compatibility conditions 
in P re i for this case. Similiarly, for j = the regularity compatibility at the boundary is 
needed if kq > 1 — [/, + 1/p, i.e., if dtp has a trace at t = 0. For the existence of a solution 
(u, p) G E ni/i x Ep jAt , it is therefore necessary that the data in (|1.5p belong to T> re \. 

(II) Let us show that (f,g,uo,po) G T> re \ is also sufficient for the existence of a unique 
solution (n, p) G E MiAt x E Pj/J of (|1.5|) . Uniqueness follows as in Step I of the proof of Lemma 
13.41 For the existence of the solution (u,p), note that it suffices to consider small T > by 
a standard compactness argument. For simplicity, we assume that is bounded. The case of 
unbounded £1 requires minor modifications. 

We cover f2 by a finite number of open balls Bi such that Bi n T = for i = 1, ...,Np and 
Bi n r 7^ for i = Nh + 1, Np, where Np, Nh G N. We further take a smooth partition of 
unity ifii for subordinate to this cover. Let (u, p) G E UjA1 x E P)M . Now, (u, p) solves (|l,5p if 
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and only if (m,pi) = (ipiU,ipip) satisfies 

d t Ui + Aui = ipif + [A, tpi]u, in f2 n B i: t G J, 

dtPi + B Ui + Copi = 7p i go + {Bo,ipi]u + [Co,'ipi}p, on T n Bi, t G J, 

BjUi + CjPi = ipigj + [B ,ipi]u+ [C ,if)i]p, on mSj, t G J, j = l,...,m, (4.1) 

u»|t=o = ipiU , in n Bi, 

Pi\t=o = i>iPo, on r n Bi, 

for all i = 1, ...,Njj, where [.A, ^]it = A(ipiu)—ipiAu. For i = 1, A^ no boundary conditions 
are involved in (|4.ip . We extend the coefficients of A outside Bi to R n such that (SD) is still 
valid, and denote the operator with extended coefficients by A 1 . Then Ui solves (|4.ip for 
i = 1, Np if and only if it solves 

d t Ui + A { Ui = ipif + [A, ipi]u, in R n , t G J, 

t*i|t=o = V^o, in M n . (4.2) 

Due to the continuity of the top order coefficients of A, the top order part of the operator A 1 
is a small perturbation of a homogeneous constant coefficient operator satisfying (E) if the ex- 
tension of the coefficients is appropriate, provided T and the radius of Bi are sufficiently small. 
Poincare's inequality in the L PjAt -spaces (Lemma 2.12 in [16J) allows to estimate lower order 
terms with constants decreasing to as T — > 0, see Lemma 1.3.13 of [14J. Using Lemma f3. 11 we 
can now solve f|4. 2() by a straightforward fixed point argument. We thus obtain a continuous 
solution operator C F : E 0jjU (J x R") x X u ^(R n ) -»• E M;M (J x R n ) for (H2J. It follows that 

Ui = C l F (ipif + [A, Tpi]u, ipiUo), i = 1, ...,N F . 

Observe that the commutator terms are of lower order. For i = Np + 1, ...,Njj the boundary 
conditions in (|4.ip are present. We choose the Bi so small that we have a chart <pi for T 
with domain Bi associated to some Xj G T. Denoting by $j the corresponding push- forward 
operator, i.e., <&iV = v o ip^ 1 , we obtain that (ui,pi) solves (|4.ip if and only if (vi,<Ji) = 
(®iUi, $ipi) solves 

dm + {<&iA<&~ l )vi = $>i(7Pif + [A,ipi]u), in R™ n ipiiBi), 

dtcii + (Siflofcr 1 )^ + Cfai = $j(^ffo + [B ,^i]u + [C ,^i]p), on R™" 1 n Pi(-B<), 

{QiB&i^Vi + Cf<Ti = ty i9j + [B is + [Cj, i>i\p) , on R™" 1 n ^(^), 

Ui|t=o = ^iV^o, in R+ n ifi(Bi), 

ai\ t=0 = ^i^ip , on R™" 1 n ^i(Si), 

for t G J and 7 = 1, ...,m. Recall that Cf 4 denotes the local representation of Cj with respect 

j J 

to the coordinates g i corresponding to tp^. According to Theorem 10.3 of |23| . at t G J and x% 
the principal parts of the operators <&iAQ~ and &iBj&~ are given by 

A^^D), B n (t,Xi,0 T v{Xi) D), 

respectively. Extending now the coefficients of the transformed operators QiAQ^ 1 , QiBjQ^ 1 
and Cj 1 such that (SD), (SB) and (SC) remain valid, we obtain that (<&iUi, &iPi) solves a half- 
space problem with operators that are either of lower order or small perturbations of constant 
coefficient operators satisfying the conditions of Lemma 13.41 As for the full-space case, if T 
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and Bi are sufficiently small, then a continuous solution operator C l H exists for this half-space 
problem, which maps the relevant data space continuously into E„ „(J x WV) x E PiAt (J x W} ). 
For i = Np + 1, A# we thus obtain 

(v,i,pi) = §J x RiC H (§i(Tl)if +[A,i/}i}u),^i(ipig+[B,ipi]u+ [C,i/)i]p),<$>i>ipiUo,$iil)ipo) . 

Here Ri is the restriction to Wl (~l ipi(Bi) and we have set B = (Bq, ...,B m ) and C = (Co, C m ) 
for simplicity. Again, the commutator terms are of lower order. 

(III) We next choose smooth functions (pi, i = 1,...,Njj, satisfying <pi = 1 on suppV'i and 
supp^j C Bi. The above considerations show that if (u, p) solves (|1.5p then it is a fixed point 
of the map Gf t g tU0 ,p ( u iP) := Yli^i^ii Pi) on the complete metric space 

Z U0 , P0 ■= {(u,p) G E„ iM x E P)At : u| t =o = «o, p|i=o = Po}- 

We remark that Z UOj p ^ s nonempty by Lemma 4.4 of |16| and Lemma 3.2.2 of [13]. Since the 
operators in the arguments of C l F and C l H are of lower order, one can show that for all data 
(f*,g*,u*,p*) G P rel the map G f* ,g* ,uq,pq has indeed a unique fixed point on Z u *^p*, making 
T and Bi once more smaller if necessary. Another fixed point argument yields for given data 
(f,g,UQ,po) € P re i the appropriate auxiliary data (/*, g* , Uq, Pq) € T> Te \ such that the fixed 
point of Gf*,g*,u^,p^ is the solution of (|1.5p . 

(IV) To finish the proof, note that the continuity of the resulting solution operator C re \ for 
(|1.5p is a consequence of the open mapping theorem. Moreover, the norm of £ re i restricted 
to "D° el is uniform in T due to an extension argument. It uses the extension operator from 
Lemma 2.5 of |16| for the oW^-spaces over J to the half-line, whose norm is independent of 
the length of J. ■ 
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